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ABSTRACT 

Ch ■ 

In this paper we propose a new strategy for gravitational waves detection from coalescing binaries, using IIR 
Adaptive Line Enhancer (ALE) filters. This strategy is a classical hierarchical strategy in which the ALE filters 

■ have the role of triggers, used to select data chunks which may contain gravitational events, to be further analyzed 
I | with more refined optimal techniques, like the the classical Matched Filter Technique. After a direct comparison 
of the performances of ALE filters with the Wiener-Komolgoroff optimum filters (matched filters), necessary 
to discuss their performance and to evaluate the statistical limitation in their use as triggers, we performed a 

• series of tests, demonstrating that these filters are quite promising both for the relatively small computational 
^vq ° power needed and for the robustness of the algorithms used. The performed tests have shown a weak point of 
ALE filters, that we fixed by introducing a further strategy, based on a dynamic bank of ALE filters, running 
simultaneously, but started after fixed delay times. The results of this global trigger strategy seems to be very 
promising, and can be already used in the present interferometers, since it has the great advantage of requiring a 
quite small computational power and can easily run in real-time, in parallel with other data analysis algorithms. 
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£jV 1. INTRODUCTION 

Gravitational Wave (hereafter GW) detection is certainly one of the most challenging goals for today physics: a 
? *~j ■ very strong proof in favor of the Einstein General Relativity description of phenomena related to the dynamics 
of gravitation and the opening of a completely new channel of information on astrophysical objects. 1-3 Ground- 
based and space detectors are now operational or will be operational in the next years, with different bands and 
sensitivities. The VIRGO/LIGO 4 ' 5 network of ground-based kilometer-scale laser interferometer gravitational 
wave detectors will probably be the key to open up a new astronomical channel of information in the frequency 
band 10 Hz to 10 kHz. In addition, when the proposed 5 million kilometers long space based interferometer LISA 6 
flies, another window will be opened in the frequency band 10~ 5 -r- 1 Hz. In particular, the sensitivities of the 
operational interferometers are becoming very interesting, producing reasonable hopes in the researchers that the 
first detection of gravitational waves is becoming very close. Therefore, the efforts in studying and developing data 
analysis methodologies and techniques are increasing, such as the tests performed on real data (interferometers 
outputs), allowing the researchers testing not only the optimality but also robustness and feasibility of such 
techniques. Within this research context, we focused our attention on GW signal detection from coalescing 
compact binaries (neutron stars (NS-NS), black holes (BH-BH) or a mixed (NS-BH)), very promising transient 
GW sources, which may also allow precise measurements of the masses of the objects, of the spins and, in the 
case of neutron stars, of the radii, too. 7 But although the interferometers seem to be sensitive enough to detect 
this class of sources, nevertheless the problem of GW signal analysis has not yet globally solved. For this task 
many solutions are being evaluated, in particular for what concerns the choice of suitable data analysis techniques 
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taking into full account the shape of the expected signal, the noise of the detector and the available computing 
power. 8 

Many efforts have been made for the development of special data analysis techniques to enhance the signal-to- 
noise ratio (SNR) of the expected GW signals. Presently, the most credited algorithm is the well known matched- 
filtering technique, based on the correlations of the detector output with templates of the expected signal. 9 ' 10 
But, although very simple in principle, its practical application requires an exact theoretical knowledge of the 
shape of the expected signal, function of the unknown parameters describing the coalescing binary and the 
evaluation of the correlation of the detector output with several thousands of templates, requirements very 
difficult to satisfy for coalescing binary GW signals. 11, 12 

The shape of the GW signal can be obtained by computing the gravitational radiation field generated by 
a system of two point-masses moving on a practically circular orbit. The solution of this problem requires the 
calculation of the gravitational radiation field to a very high order in terms of a post-Newtonian expansions, 
being coalescing compact binaries strongly dominated by relativistic effects. The large number of templates 
necessary for data analysis using matched-filtering technique arises problems due to the great computing power 
needed to perform this task on-line. 8, 11 In fact, as a consequence of the large band of interferometric detectors 
(some kHz), sampling rates of the order of 20 kHz are foreseen, resulting in a huge amount of data/day to be 
analyzed on-line (~ 10 GByte/day). Moreover, the computational cost depends on the number of parameters 
used for the phase approximation, on the accuracy of the sampling of the likelihood function (in connection with 
the capability of recovering weak signals) and on the chosen frequency band. For example, a direct application 
of the matched-filtering technique to the VIRGO antenna requires a computing power of at least « 0.3 Teraflops 
starting from a minimum mass of 0.25M© with a SNR recovery of 90%, including the computing power necessary 
for the production of the templates. 13, 14 These estimates confirm that already at the stage of the first post- 
Newtonian correction large parallel computers are necessary for an on-line analysis of these signals. Of course, 
the analysis of such a large amount of information could be made off-line, but it is clear that great advantages 
may derive from a preliminary selection of the on-line data frames which may contain a GW signal, not last the 
fact that the selected data chunks will be object of more refined analysis, including also cross-correlations with 
data taken by other running interferometers. 

Starting from the estimate of the huge computational cost of the data analysis for VIRGO, we decided to 
explore alternative approaches to try to solve this problem that can be used also in parallel to the existing 
ones. For this reason, also taking into account the relevance of a GW first detection, we decided to implement 
a hierarchical technique, based on fast on-line sub-optimal filters to be used as triggers (first level), in order to 
extract interesting data chunks, followed by an optimal technique like the matched filters one. 15 Our efforts 
were therefore aimed to the implementation of a rough analysis with the following characteristics: minimum 
signal losses with respect to the one using the matched Wiener-Kolmogorof filter, robustness against false alarm 
detection and, last but not least, need of low computational powers and possibility of running in real-time. 

2. IIR ADAPTIVE LINE ENHANCER 

In a previous paper, 16 we tested an adaptive filter of the class of IIR (Infinite Impulse Response) Adaptive 
Line Enhancers (IIR ALE). The ALE algorithm is a technique designed to approximate the SNR gain obtained 
by the matched filter solution for a sinusoidal signal. The advantage of using ALE filters is that no a priori 
knowledge of the signal parameters (sinusoidal frequencies, amplitudes, phases or even the number of narrow- 
band components) is required. Therefore, center frequency band-pass filters can be used when the signal to detect 
is narrow-band and buried in a wide-band noise. And a slow frequency sweep very well approximate this case. 
In fact, instead of searching for the optimum Wiener Filter, that requires the knowledge of a general frequency 
response to be adapted to the signal, the basic idea of this procedure is that of starting from a band-pass response 
and then optimize its parameters. Something similar is done in Matched Filter technique, where the occurrence 
of a signal must be detected, in the hypothesis of known signal waveform, but with unknown parameters. 

The optimization procedure may be carried out in the following way. Let us suppose that the observed 
signal is x n — s n + w n , where s n is a band-limited signal (eventually a single sinusoid) and w n is a white noise 



uncorrelated with s n . Let Y(z) be the output of a digital filter with transfer function H(z) in the Z-domain: 

Y(z) = H{z)X(z) = H(z)[S(z) + W(z)} (1) 
The error after filtering in the time domain will then be 

dn Vn S n (2) 

and in the Z-domain 

D(z) = Y(z) - S(z) = H(z)W(z) + [H(z) - l]S(z). (3) 

If signal and noise are assumed to be uncorrelated, referring to the above formula, the power spectrum of the 
error signal is 

S d (w) = \H(e? u ) - if S s (uj) + \H(e^)\ 2 S w (u). (4) 
The rms error can then be evaluated by inverting this formula, leading to 

E R} = ±J \H{en - if S s {u)du + 

i- J \H(en\ 2 S w (u)duo. (5) 

Then choosing the root mean square error as the cost function to be minimized, the optimization will lead to 
the optimal filter transfer function H(e^). Supposing that our simple model, consisting in a sinusoid immerse 
in gaussian noise, can be described in the frequency domain by 

S s {e ju ) ^2tt5{lj-u ) and S w (u) * o*, (6) 

then the root mean square error can be written as 

^K} = ^l^ )-l| 2 + ^ / \H{en\ 2 du (7) 

which is minimized by the two conditions 

J \H{e^)\ 2 dw = min! with H(e jw °) = 1 (8) 

2tt 

thus leading to a sharp band-pass filter. 

On the basis of the above considerations, the idea embedded in a band-pass Adaptive Line Enhancer (ALE) 
is that if H(e JU} ) is a parametric band-pass transfer function, described by one or more parameters, then it is 
possible to find the best value of these parameters by calculating the minimum of the above error. If the analytical 
structure of the error function prevents from a closed form solution, we may rely on a numerical optimization 
both in the time and in the frequency domain. The above procedure is well suited to the ideal case, when all the 
parameters are well known. In the real case the signal to be detected has unknown or time varying parameters. 
Therefore, an adaptive implementation of the filter becomes necessary to detect the signal and to trace its time 
varying features. 

2.1. Optimization of ALE Parameters 

Suppose that the signal model is the sine wave plus (white) noise model 



x n = a sin(oj n + 0) + w. 



(9) 



input of an adaptive filter with a band-pass parametric structure. The optimization must be then performed 
looking for a maximum of its output power a 2 .. The performance function in this case is 



E[\y n \ 2 } = Z- 1 (S y (z))\ n=0 (10) 

Then since 

Z- 1 (S v {z)) = ^-js v {z)z^dz (11) 

then in the case of a linear filter with frequency domain response H{e^ ul ), the expected power of the filter output 
can be written as 

E {\Vn\ 2 } = ^ J S x {w)\H{e^)\ 2 dw (12) 

In the simple case of sine wave or band-pass signal 

S a (u;) S ±=8(u> - uj ) (13) 

within uncorrelated white noise, S w (u>) = <r 2 , we have 



E 



{\y n \ 2 } - y \H(e^)\ 2 + || £ \H(ju;)\ 2 (14) 



It is necessary to evaluate the performance function, in order to study the features of the performance surface 
of the system, and to define the most suitable algorithm to reach the maximum. An unimodal surface will 
obviously simplify this search, avoiding the need for the use of a global search algorithm. Many implementations 
are reported in literature, characterized by different structures. Within this context two are the most relevant 
parameters: convergence speed and tracking capabilities, that are clearly contrasting requirements. An improve- 
ment may be achieved if the parameters of the optimization include both filter center band and the bandwidth, 
at least in the final stage of the convergence. This may improve the main feature of the method which is the 
capability to provide a running search of the center frequency, improving the detection of the instantaneous 
frequency. 

2.2. Choice of the prototype band-pass filter 

In principle any adaptive band-pass filter is suitable for our purpose. Anyway it may be convenient to properly 
shape the band in order to exactly match the shape of the signal in the frequency domain, since it will act as a 
weighting function in the frequency domain. Therefore, in this case it is convenient to use a filter prototype with 
the selected shape in the frequency domain, to which apply a parametric frequency transformation, in order to 
obtain the final band-pass filter. 

We implemented the second-order band-pass Butterworth ALE filter 19 as 



H(z) = (l-r 2 ) 



W t z/{r + r 2 )-l 



(15) 



where W t — 2rcos(27r/ ) is the only parameter dependent on the center-frequency / , while r is a fixed design 
parameter related to the filter bandwidth B by the relation B = F s (1 — r)/2, where F s is the sampling frequency. 
The recursive algorithm for the filter can be expressed by: 

yt = Wt(^^jx t ^-{l-r 2 )x^ 2 +W t y t ^-r 2 y t ^ (16) 

«t = = (~~T~) Xt ~ x + w t a t~i + yt-i - r 2 a t -2 (17) 

Rt+i - iyRt+a 2 (18) 
W t+1 = Wt + ^p^ (19) 



where v G (0, 1) is the forgetting factor, introduced for the recursive computation of the normalizing factor R, 
a t is the instantaneous gradient and \i a is the scalar adaptation step size. Stability monitoring of these filters is 
very simple, being performed in order to verify the condition \Wt\ < 2r. 

Let us assume a coalescing binaries signal (chirp) injected in AWGN (Additive White Gaussian Noise) with 
a theoretical SNR% n * . The output of ALE filter has an improvement of 

SNR\ LE \ _J\H(z)\^ ss (dz/z) (2Q) 



SNRl J ~ 2nja*Bf 

where $ ss is the power spectral density of the signal, <7„ is the input noise power and Bf is the equivalent noise 
bandwidth. When the ALE filter tracks the input signal, / tends to coincide with the signal instantaneous 
frequency and, therefore, the equation 20 can be written as 15 

SNR ALE = SNRl ■ (21) 

From equation 21 we can state that, if the condition (1 — r) « 1 is satisfied and ALE filter closely tracks the 
signal, there is an improvement of SNRf n by a factor of (1 + r)/(2(l — r)) simply because ALE is a band-pass 
filter with a adaptive center frequency 

2.3. ALE thresholds: probability of false alarm and probability of detection 

Using an ALE filter for detection, a very delicate problem is the choice of the threshold T, to discriminate 
between the alternative hypotheses of absence or presence of a signal at an assigned level of confidence. Since 
we use an ALE filter, we can assume that we have an approximation of the instantaneous signal amplitude, 
while the signal phase is unknown. In these hypotheses, in presence of only gaussian noise with variance cr^, the 
output power of the ALE filter is ajj — 2er^(l — r)/(l + r), and it can be modeled with a random variable having 
a Rayleigh distribution with spread parameter c^ 9 : 

Po(x) = -J- cxp j-^r} (22) 

a N I N ) 

If a signal is present and is tracked, then the output power statistic becomes a Rice distribution centered on the 
signal energy £ , 9 and with spread parameter cr| = 2(7^(1 — r)/(l + r) = crj^: 



x T / xE 



x 2 +£ 2 



p 1 (x) = —I [ — )e 2a s with.x>0 (23) 

(7 a \ (7 . 



where Jo(-) is the modified th order Bessel function of first kind. The false alarm probability, assuming inde- 
pendent samples, can be written as 



Pfa = J p (x)dx = exp j- 



2< 



(24) 



that is a suitable closed form for the evaluation of the optimum detection threshold T that can be calculated as 



T = 1 /-2a2 r logP /Q = <7 M /-4 ( — ) logP /a (25) 



Concerning the probability of detection P4, assuming that the signal is present and tracked, then we have 

r+00 

P d = J t p\(x)dx — Q (SNRale,T/ SNRale) (26) 

"The definition of signal to noise ratio we used is SNR% n = J 2 % where E = ^2 i xf(t) is the signal energy, a 2 is the 
noise variance and F s is the sampling rate. 



where Q is usually referred to as the Marcum Q function 



Q(a,b) 



a -\-x 

xe 5 I (ax) dx 



(27) 



We made several synthetic experiments to test the proposed algorithm. We used the waveform generated by 
coalescing binary stars with masses 1.4 — 1.4M Q at the 2.5 post newtonian expansion. The noise source has a 
Gaussian probability distribution in the frequency band 50 -j- 5000 H z. The ALE filter has a starting frequency 
of 50 Hz and a fixed band of 25 Hz corresponding to r = 0.990. The empirical ROC curve of the ALE filter, 
necessary to characterize the performance as a sub-optimal detector, is shown in figure 1. It is immediately clear 
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Figure 1. The empirical ROC curve of single ALE filter, evaluated at different SNRale, is shown. 

that a high probability of detection within a high false alarm probability is obtained also at very low signal- 
to-noise ratio, demonstrating, in this way, the feasibility of using an IIR ALE filter as trigger with an assigned 
confidence level. In fact, choosing a higher false alarm probability, we may use an ALE IIR like an efficient 
trigger also at very low SNR. In figure 2 an example of ALE filter output is shown: from the top it can be seen 
a data set of binary signal injected in gaussian noise with SNRale — 10, the binary signal, and the frequency 
estimated by ALE filter. 

2.4. Parallel implementations 

An accurate analysis of the ALE performances shows that in the low SNR case only the choice of a starting center 
frequency of the ALE filter very close to the frequency to detect may allow taking advantage of the gradient 
procedure thus reducing the effect of high level noise. In our case, although the transition of the gravitational 
wave chirp through all the frequencies of the sensitivity band of the interfcromctric detectors is theoretically 
predictable, the time at which the gravitational wave chirp passes through the starting center frequency of the 
ALE filter is not predictable. Actually this may be a problem for the ALE filter. In fact, when a GW chirp 
reaches the design starting center frequency of the ALE filter, if its present center frequency has moved from its 
design one due to the input noise effects, then the probability of detection largely decreases. 

In order to keep this probability the highest possible, closest to the maximum theoretical one, it is important 
to synchronize the design starting center frequency of the ALE filter with the GW signal frequency. A possible 
solution to this problem is to use parallel techniques, like the one of building a bank of ALE (band-pass) filters, 
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Figure 2. An example of Ale filter output is shown. From the top, a data set of binary signal injected in gaussian noise 
with SNRale — 10, the binary signal and the frequency estimated by ALE filter, are shown. 



with a common single design starting center frequency, but starting at different times, sequentially and with 
fixed delay times. The number of ALE filters of the bank is defined at the beginning of the procedure. When 
the bank is filled, and a new ALE filter is started, then the oldest one is eliminated. 

The basic idea on which this procedure is based on, is clearly that of overcoming the lack of synchronization 
among the ALE filters at start-up and the possible presence of a GW signal. Of course, the shortest is the 
delay time, the best is the synchronization. In the limit it is possible to make this delay time coincident and 
synchronous with the inverse of the sampling frequency of the data acquisition system. In fact, due to the low 
computational power required it is possible also to think of a bank of millions of ALE filters, running in parallel 
for a real-time data analysis. This solution allows a faster identification of the frequency to be detected, coupled 
with a more robust signal identification and tracing. In fact, if the ALE filters are very close each other, then in 
presence of a GW signal it is obvious that a quite large number of them should detect the signal, allowing also 
to perform a better detection. What need to be still to be completed is the development of a real-time decision 
unit, necessary to statistically interpret the data from the parallel filter bank and to control the center filter 
frequencies. 

All the numerical experiments carried on by us, in low SNR conditions, exploited this parallel approach, and 
produced very good results. The empirical ROC curve of the parallel ALE filter, necessary to characterize the 
performance as a bank of sub-optimal detectors, is shown in figure 3. As it can be seen, an improvement of the 
performances is archived, respect to the single ALE filter. 

Finally, for completeness, a thoroughly parallel structure could also be used in order to trace the error surface, 
for example by sweeping over the entire frequency range, thus avoiding possible looping on local minima, and 
then providing robustness to the detection. 

3. CONCLUSIONS AND FUTURE IMPROVEMENTS 



We have analyzed and discussed the performances of an IIR ALE filter to perform real-time triggering on the 
output data from interferometric antenna, like the VIRGO one, showing that it is suitable for GW data analysis 
from coalescing binaries. This filter is based on adaptive modification of a band-pass filter parameters, changed 
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Figure 3. The empirical ROC curve of Parallel ALE filter, evaluated at different SNRale, is shown. 

according to a suitable cost function, within the hypothesis that the signal to detect is narrow-band (even if the 
center band is changing along the time) buried in a wide-band noise. For this task we performed some choices, 
analyzing the performances of 2 nd order direct implementation of IIR adaptive filters, either in the case of a 
single filter or in the case of a parallel bank of filters. From the analysis of the detection probability, we obtained 
the results that are shown in figure 1 and 3. 

The simulations we performed demonstrate that on-line triggering on gravitational wave interferometric 
experiments is feasible without a huge or complex computing power and without the necessary a priori knowledge 
of the signal waveform. Coupling a parallel process of different algorithms will surely provide a better and robust 
on-line tagging of binary coalescence signals, although further studies arc necessary to test in greater detail the 
parallel version of ALE triggers. 
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